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Abstract 

In this paper we have investigated an LRS Bianchi I anisotropic cosmo- 
logical model of the universe by taking time varying G and A in the pres- 
ence of bulk viscous fluid source described by full causal non-equilibrium 
thermodynamics. We obtain a cosmological constant as a decreasing func- 
tion of time and for m, n > 0, the value of cosmological "constant" for 
this model is found to be small and positive which is supported by the 
results from recent supernovae observations. 
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1 Introduction 



The conventional theory of evolution of the universe includes a number of dissi- 
pative processes. Dissipative thermodynamics processes in cosmology originat- 
ing from a bulk viscosity are believed to play an important role in the dynamics 
and evolution of the universe. Misner[T] suggested that large-scale isotropy of 
the universe observed at the present epoch is due to action of neutrino viscosity 
which was not negligible when the universe was less than a second old. A number 
of processes responsible for producing bulk viscosity in the very early universe 
are such as the interaction between radiation and matter :2; , gravitational string 
production (3J 0] , viscosity due to quark and gluon plasma, dark matter or par- 
ticle creation jS[ It is important that each dissipative process is subject to as 
detailed analysis as possible. However, it is also important to develop a robust 
model of dissipative cosmological processes in general, so that one can analyze 
the overall dynamics of dissipative without getting lost in the details of partic- 
ular complex processes. In requirements of such a model MaartensJT] pointed 
out that the model should (i) be causal and stable, and (ii) provide a constant 
relativistic thermodynamics in the 'conventional' post-inflationary regime. 

In order to study these phenomena, the theory of dissipative was first de- 
veloped by EckartjHj and subsequently an essential equivalent formulation was 
given by Landau and Lifshitz|5|. But Eckart theory has several drawbacks in- 
cluding violation of causality and instabilities of equilibrium states. Readers 
interested in the general theory of causal thermodynamics are urged to consult 
the excellent survey report of Maartens|7| and Zimdahl and references cited 
therein. A relativistic second-order theory was found by Israelii] and devel- 
oped by Israel and Stewart [T2*]. The advantages of the causal theory are as 
follows [T5|: (1) For stable fluid configurations, the dissipative signals propagate 
causally. (2) Unlike Eckart-type's theory, there is no generic short- wavelength 
secular instability in causal theory. (3) Even for rotating fluids, the pertur- 
bations have a well-posed initial value problem. Therefore, the best currently 
available theory for analyzing dissipative processes in the universe is the full 
(i.e. non-truncated) Israel-Stewart causal thermodynamics, which we consider 
in this work. 



In recent years, models with relic cosmological constant A have drawn con- 
siderable attention among researchers for various aspects such as the age prob- 
lem, classical tests, observational constraints on A, structure formation and 
gravitational lenses have been discussed in the literature (see Refs.|14 | — |16 p . 
Lindey^Z] has suggested that cosmological "constant" may be considered as 
a function of temperature and related to the spontaneous symmetry breaking 
process. Therefore, A should be a function of time in a homogeneous uni- 
verse as temperature varies with time. Some of the recent discussions on the 
cosmological constant "problem" and on cosmology with a time- varying cosmo- 
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logical constant are given by Ratra and Peebles JS|> Dolgov[T^| — |2T]. Sahni and 
Starobinsky 22 , Peebles[2U, Padmanabhan pi| and Vishwakarma|25| . Recent 
cosmological observations suggest the existence of a positive cosmological con- 
stant A with the magnitude A(Gh/c 3 ) « 10~ 123 . There are several aspects of the 
cosmological constant both from cosmological and field theoretical perspectives. 
Presently, determination of A has become one of the main issues of modern cos- 
mology as it provides the gravity vacuum state and make possible to understand 
the mechanism which led the early universe to the large scale structures and 
to predict the fate of the whole universe. The cosmological "constant" can be 
measured by observing quasars whose light gets distorted by gravity of galaxies 
that lies between the quasars and Earth. Krauss and Turner pB] have mentioned 
that as A term dominates the energy density of the universe, cosmologists are 
correct in their attempt to evoke it once again for better understanding of both 
the universe and fundamental physics. 

In the last few decades there have been numerous modifications of general 
relativity in which gravitational "constant" (G) varies with time|27]. Consider- 
ing the principle of absolute quark confinement, Der Sarkissianj^HI has suggested 
that gravitational and cosmological "constant" may be considered as functions of 
time parameter in Einstein's theory of relativity. A number of authors 29 — 40 
have considered time- varying G and A within the framework of general relativ- 
ity. 

Motivated by the fact that bulk viscosity, gravitational and cosmological 
"constants" , are more relevant during early stages of the universe, in this pa- 
per, we have considered the evolution of a LRS Bianchi I model with bulk viscous 
fluid in full causal non-equilibrium thermodynamics, in presence of time- varying 
gravitational and cosmological "constants" . 



2 THE Basic Equations 

A locally rotationally symmetric (LRS) Bianchi I space-time may be represented 
by the line element 

ds 2 = dt 2 - A 2 dx 2 - B 2 {dy 2 + dz 2 ), (1) 

where metric potentials A and B are depending on cosmic time t only. 
The Einstein's field equations with variable G and A are given by 

Rij - \R9H + A (%y - -87rG(t)T y , (2) 

where Rij is the Ricci tensor; B = g 1 -* Rij is the Ricci scalar; and Ty is the 
energy-momentum tensor of cosmic fluid in the presence of bulk viscosity given 
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by 



Tij = (p+p + n)uiUj - pgtj. 



(3) 



Here p, p and II are the energy density, equilibrium pressure and bulk viscous 
pressure respectively and u 1 is the flow vector satisfying the relations u l Ui = 1. 
The Einstein's field equations for the line element Q lead to the following 
set of equations 

|j = 87rG(t)p + A(t), 



2AB 
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B 



AB 
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j = -8nG(t)[p + U}+A(t), 



If + = -87rG(t)[p + II]+A(i). 
A combination of equations Q - JBJ yield the continuity equation 



A = 0. 
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From equations Q and ©, we get 

A = SirGp. 



(8) 



(9) 



For the full causal non-equilibrium thermodynamics, the causal evolution equa- 
tion for bulk viscosity is given by|7] 



rn + n 




r n 




(10) 



where T > is the temperature, £ the bulk viscous coefficient and r > the 
relaxation coefficient for the transient bulk viscous effect (relaxation time i.e. 
the time which system takes in going back to equilibrium once the divergence 
of the four velocity has been switched off). For r = 0, equation IjlOjl gives evo- 
lution equation for the non-causal theory. For e = 0, we get causal evolution 
equation for truncated theory, which implies a drastic condition on the temper- 
ature, while e takes value unity in full causal theory. 
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3 The Model 

Equations (0 and JBJl yield 

B B 2 A AB 

B + B^-A-AB-°- <M 

In order to find exact solutions of the field equations, we require more physically 
plausible relations amongst the variables. Considering a power law relation 
between A and B as A ~ B n , ea. ljTTJl suggests 

(12) 

a') 

where Ao and Bq are the values of A and -B at present time t — to and n is a 
positive constant. 

By use of l|T^|l and (fH3|) . equations ® - ®, reduce to 

2T7 1 

8ttGp + A, (14) 



(n + 2) 2 t 2 
2n + l 1 



(n + 2)2t2 ^(p + ID-A. (15) 

Further following 2 Tj — US], we assume G = t m . Hence equations (0 and 11411 
suggest 

: m . 2mn 1 . . 

A A = -7 ^fo^) ( 16 ) 

which is a linear equation and it has solution 

A — 2mn 1 

~ (n + 2) 2 (m + 2)i 2 ' ( ' 

From equations l|14|l and l|17(l one can easily obtain expression for energy density 
in terms of cosmic time t as 

(18) 



r 27r(TO + 2)(n + 2) 2 i m + 2 ' 

Considering the usual barotropic equation of state relating the perfect fluid 
pressure p to the energy density as 

p=( 1 -l)p, (19) 

where 7(1 < 7 < 2) is a constant, equations i|13|) and (|14|l lead to 

n = rijm + l) 1_ 

27r(m + 2)(n + 2) 2 t 2 + m ' 1 ; 
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Now, we consider |7I HU | |4"l" ] — 144 | following phenomenological widely accepted 
relations 

£ - S, p a and r=- (21) 
P 

for the bulk viscosity coefficient £ and mass density p and also for bulk viscosity 
coefficient and the relaxation time r, respectively, where £o > and a are 
constants. If a = 1, ea. (|21|l may correspond to a radiative fluid, whereas a = § 
may correspond to a string-dominated universe. However, more realistic models 
are based upon a in the region < a < |. 



Using (|21|l . equation UlUfl on integration yields 



T = T exp 



- e 9(t) 



exp 



117 R 3 



(22) 



where f(t) and g(t) are anti-derivatives of i and respectively. 

With the help of equations 118(1 . (|20|l and J2J, we obtain the expressions for 

f(t) and g(i) as 

to 



where 



/(*) = 
9(t) 

1 



'.+2) (a 
1 



m(m + l)t T 



1 - (m + 2)(a-l) 



27r(m + 2)(n + 2) 2 



(23) 
(24) 

(25) 



We observe from eq. (|17() that the cosmological constant in this model is de- 
creasing function of time and it approach a small value as time increases ( i.e., 
the present epoch). For m, n > the value of cosmological "constant" for this 
model is found to be small and positive which is supported by the results from 
recent supernovae observations recently obtained by the High - z Supernova 
Team and Supernova Cosmological Project ( Garnavich et aZ.;0F] Perlmutter et 
a/.;gnj Riess et a/.; 03 Schmidt et al. |H]). 

Using the observational values G/G = 10 _11 yr _1 and H p = 7.5 x 10~ 11 yr~ 1 , 
and the relation of the present age of the universe with Hubble parameter 
(H p t p ~ |), relation G — t m suggest 



22.5 



This clearly shows that gravitational parameter G turns out to be an increasing 
function of time. 

From eq (|18fl . it is observed that the energy density p is decreasing with evolution 
of the universe. For all m > —2, n > 0, the energy conditions are satisfied. 
It is also observed from eq (|2lH) that bulk viscous pressure n decreases with 
time. When t -> 0, n -> oo. When t -> oo, n -> 0. 
In this model, expressions for expansion and shear are 



= — and 



a 2 = 



2(n-l) 2 
3(n + 2) 2 t 2 ' 



(26) 



G 



^ = 47r(n-l) 2 (m + 2) fm 
p 3n 

Equation (|27[) clearly shows the effect of time varying G on the relative anisotropic 
For, —2 < m < 0, relative anisotropy is decreasing. Further, it can be observed 
that cr 2 oc 9 2 , which indicate that the model does not approach isotropy for 
large value of t. 
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